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两 类 伽 马 暴 X 射线 余辉 光 变 特性 的 比较 研究 
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摘要 : 伽 马 射线 暴 ( Gamma-Ray Bursts, GRBs) 的 分 类 是 伽 马 暴 领 域 未 解决 的 热点 问题 ， 
为 了 解决 这 个 问题 ， 我 们 搜集 了 Swift 观测 具有 缓慢 衰减 和 随后 衰减 特征 的 又 射线 余辉 光 变 
曲线 样本 (180 个 已 知 红 移 的 伽 马 暴 )， 其 中 长 伽 马 暴 (Long Gamma-Ray Bursts ，LGRBs ) 172 
个 ， 短 伽 马 暴 (Short Gamma-Ray Bursts, SGRBs)8 个 ， 并 选择 平滑 摇 折 需 律 函数 使 用 马尔 可 
夫 链 蒙特 卡 罗 ( Markov Chain Monte Carlo, MCMC) 方法 拟 合 了 这 些 X 射线 余辉 光 变 曲线 ， 获 
得 光 变 特性 参数 后 ， 分 析 两 类 暴 的 光 变 特性 参数 a (ARERR), а, (ЕЖА). 
也 (静止 系 下 缓慢 衰减 结束 时 间 ) 的 分 布 是 否 存 在 差异 以 及 两 类 暴 是 否 遵循 相同 的 相关 关系 。 
研究 发 现 ，(]1) 两 类 暴 的 光 变 特性 参数 ai о, 和 也, 的 分 布 范围 、 中 值 和 平均 值 都 非常 接近 ， 
没有 显示 出 明显 的 不 同 ， 随 后 的 K-S 检验 结果 进一步 表明 ， 这 两 类 暴 的 光 变 特性 参数 分 布 
没有 显著 差异 。(2) 我 们 选取 的 两 类 暴 均 遵循 如 下 相关 关系 : L.-T, LE ios ST f ЕЁ, io 
7%， 而 且 它 们 在 散 点 图 上 不 可 区 分 。 这 些 统计 结果 表明 ， 两 类 暴 的 这 些 光 变 特 性 是 相似 的 ， 
因而 ， 我 们 推断 两 类 伽 马 暴 具 有 同一 物理 起 源 ， 它 们 的 喷 流 结构 也 可 能 相同 。 
关键 词 : лаж; X 射线 余辉 ; 光 变 曲线 拟 合 ; 数据 分 析 
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伽 马 射 线 暴 ( 简称 伽 马 暴 ) 是 宇宙 中 短暂 、 强 烈 的 伽 马 射线 闪光 ， 持 续 时 间 从 几 毫 秒 到 几 千 秒 不 
等 天 1。 根据 持续 时 间 和 光谱 硬度 比 的 分 布 ， 伽 马 暴 可 以 分 为 两 类 : 长 伽 马 暴 (Tu<2s， 通 常 持续 
20~30 s) 和 短 伽 马 暴 (SGRBs,， Tu<2s， 通 常 持 续 0.2~0.3s)!5 。 通 常 认为 ， 长 伽 马 暴 起 源 于 大 质量 
恒星 的 塌 缩 '“ ,得 伽 马 暴 起 源 于 两 颗 致 密 恒星 的 并 合 ， 尤 其 是 两 颗 中 子 星 (NS-NS ) O 或 一 颗 中 子 星 
和 一 个 黑洞 (NS-BH) "9 。 对 于 СЕВ 170817А 相关 的 引力 波 事件 GW 170817 的 检测 证 实 ， 至 少 有 一 
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目前 已 认识 到 伽 马 射 线 暴 包含 两 个 现象 的 辐射 阶段 : 瞬时 辐射 (瞬时 伽 马 射线 辐射 ) 和 余辉 阶段 
(瞬时 后 的 长 持续 、 宽 能 带 辐射 ) "31 。 火 球 激 波 模型 广泛 应 用 于 解释 伽 马 射线 暴 的 起 源 ! 2 。 研 究 发 
现 ， 瞬 时 辐射 由 火球 喷射 物 相 互 作用 引起 的 内 部 激 波 产生 ， 而 余辉 辐射 由 火球 喷射 物 与 星际 介质 相互 
作用 引起 的 外 部 激 波 产生 ""””。1997 年 ，BeppoSAX 卫星 凭借 精确 的 定位 能 力 首次 探测 到 伽 马 射线 
炊 的 X 射线 余辉 ， 同 时 还 检测 到 件 马 暴 的 光学 和 射电 余辉 ， 使 伽 马 暴 的 定位 得 以 精确 。 文 [26] 首 次 
测量 GRB 970505 的 红 移 为 0.835， 证 明 伽 马 射 线 暴 起 源 于 宇宙 距离 ， 后 来 甚至 观测 到 红 移 ~100” 1 。 

随 着 Swift 卫星 的 成 功 发 射 ， 研 究 人 员 发 现 了 早期 X 射线 余辉 的 许多 特殊 特征 ， 例 如 X 射线 耀 发 
和 缓慢 衰减 成 分 ”1 。 文 [28] 提 出 了 具有 5 个 阶段 典型 的 X 射线 余辉 光 变 曲线 ， 包 括 快速 衰减 成 分 、 
缓慢 误 减 成 分 、 正 常 衰 减 成 分 、 后 喷 流 拐 折 成 分 和 X 射线 炊 发 。 通 常 认 为 ,缓慢 衰减 阶段 是 一 些 额 
外 的 能 量 注入 。 这 可 能 是 由 于 长 期 的 中 央 发 动机 持续 释放 的 能 量 '””。 因 此 ， 人 们 相信 X 射线 余辉 
缓慢 衰减 可 能 有 助 于 确定 伽 马 射线 暴 不 同 中 央 引 擎 的 类 型 。 其 中 一 个 广泛 接受 的 理论 是 具有 X 射线 
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令 人 惊讶 的 是 ， 在 缓慢 衰减 阶段 结束 后 ， 人 们 还 发 现 了 另 一 种 特殊 的 缓慢 衰减 (所 谓 的 “内 部 平台 ” 
它 具 有 一 段 及 其 快速 的 衰减 (衰减 指数 通常 小 于 -3)。 一 种 公认 的 解释 是 ,“ 内 部 平台 ”是 由 于 “超大 
质量 ”中 子 星 的 持续 能 量 供应 。 中 子 星 最 终 会 在 自转 或 回落 吸 积 的 作用 下 塌 缩 成 黑洞 ， 导 致 中 央 发 动 
机 突然 关闭 > SOO) 。 与 “内 部 平台 ”相反 ， 通 常 称 为 “外 部 平台 ”的 缓慢 衰减 后 的 正常 衰减 ， 这 
可 以 很 好 地 解释 为 外 部 激 波 的 减速 。 
在 具有 X 射线 余辉 缓慢 误 减 阶段 的 伽 马 射线 暴 中 ， 研 究 人 员 发 现 了 几 个 重要 的 相关 性 。 首 先 ， 
文 [47] 分 析 了 77 个 伽 马 射 线 暴 ， 发 现在 静止 坐标 系 下 慢 误 减 阶段 的 结束 时 间 ( 了 , ,) 与 此 时 对 应 的 
射线 光度 (L,) 之 间 存 在 负 相 关 关系 ，L, x70 ”然后 ,通过 使 用 具有 已 知 红 移 的 62 个 Swift 
长 件 马 暴 样 本 ,他 们 将 其 扩展 到 余辉 与 瞬时 辐射 之 间 ， 并 发 现 L ,与 瞬时 辐射 平均 光度 Ls 也 存在 相 
关 性 ， 即 ,x 7023s。 文 [50] 进 一 步 发现 7,， 工 ,和 瞬时 辐射 的 各 向 同性 能 量 (E, ;, ) 之 间 存在 更 
紧密 的 三 参数 相关 性 。 文 [51] 建 立 了 瞬时 辐射 也 ,， 到 ,和 峰值 光度 (二 ) 之 间 的 另 一 个 三 参数 相关 性 。 
然后 由 文 [52] 证 实 扩展 了 183 个 伽 马 射线 暴 的 样本 。 文 [53] 分 析 了 50 个 在 光学 余辉 光 变 曲线 中 具有 
绥 慢 衰减 阶段 的 伽 马 射线 暴发 现 ， 样 本 具有 相似 的 三 参数 L-T-E 相关 性 。 后 来 ， 文 [54] 分 析 了 174 
рй 个 具有 缓慢 衰减 阶段 和 已 知 红 移 的 Swift MERKER, AMAT L-T-E 的 相关 性 ， 进 一 步 证 实 了 
С 这 一 相关 性 ， 并 得 出 结论 ; 1, ост. 09 p 30. 0t 
然而 ， 前 人 的 研究 并 没有 详细 分 析 和 系统 比较 具有 缓慢 下 降 阶段 两 类 暴 的 X 射线 余辉 光 变 特性 ， 
而 目前 伽 马 暴 的 分 类 问题 还 是 一 个 未 解决 的 热点 问题 。 在 这 项 工作 中 ， 我 们 从 两 类 伽 马 暴 余 辉 的 光 变 
特性 出 发 ， 系 统 收集 X 射线 余辉 光 变 曲线 数据 作为 研究 样本 ， 重 点 关注 X 射线 衰减 后 的 阶段 ,统计 
两 类 暴 平台 后 面 的 拐 折 情 况 ， 哪 类 暴 更 倾向 于 有 这 样 的 平台 。 我 们 使 用 平滑 拐 折 寡 律 函 数 模型 对 两 类 
伽 马 暴 的 X 射线 余辉 光 变 曲线 进行 拟 合 ， 得 到 缓慢 下 降 阶段 的 相关 参数 ， 从 统计 上 分 析 两 类 伽 马 暴 
的 缓慢 下 降 阶 段 光 变 特性 的 异同 ， 为 推 新 伽 马 暴 的 物理 起 源 、 暴 周 环境 、 前 身 星 、 辐 射 机 制 等 问题 提 
供 帮 助 。 本 文采 用 平坦 ACDM 宇宙 学 模型 ，H,=70.0 km +s! + Мре Ail Q, = 1-0, =0. 286 计算 D(z)。 


1 样本 选择 和 数据 分 析 


Swift 在 2004 年 成 功 升 空 后 !5 ， 很 快 就 发 现 了 伽 马 暴 X 射线 余辉 许多 有 趣 的 特性 ， 包 括 X 射线 
= Ж Ac AIZEN FEB ES, REZE Swift XRT 探测 到 有 缓慢 下 降 以 及 随后 正常 下 降 阶 段 的 伽 马 暴 。 
Б 我 们 选取 从 2005 年 3 月 到 2022 年 5 月 的 伽 马 暴 ， 所 有 伽 马 暴 都 测量 到 红 移 ， 能 量 范 围 在 0.3~ 10 keV. 

从 Swift-XRT {jl 438628 th 2S FE (https: // www.swift.ac.uk/xrt_curves/) 中 ， 我 们 根据 以 下 标准 筛选 
符合 要 求 的 伽 马 暴 : 

(1) 缓 慢 下 降 阶 段 比 较 明 显 ， 因 此 要 求 缓慢 下 降 阶 段 的 寡 律 指数 在 -1.0 到 1.0 的 范围 内 ; 

(2) 缓 慢 下 降 阶 段 的 数据 足够 丰富 ， 既 能 显示 缓慢 下 降 阶段 的 结构 ， 又 能 显示 缓慢 下 降 阶 段 的 结 
ЖЇР]; 

(3) 218 REBT EZRA TAA, [ГЕЛИ ЖЕЛ ФТ ХАИИА ШЕ, ШИ ЖИПКЕ Р A E 
发 ， 我 们 也 只 取 其 缓慢 下 降 阶段 以 及 衰减 之 后 的 数据 进行 拟 合 ; 

(4) 样 本 的 红 移 已 知 ， 这 样 我 们 就 可 以 计算 y 射线 释放 的 各 向 同性 能 量 ( 已 ,、.) ， 并 且 还 能 够 考虑 
时 间 脱 胀 效应 推导 内 在 T, 参数 。 

通过 以 上 4 个 标准 ， 我 们 最 终 获 得 了 180 个 伽 马 射线 暴 样本 ， 其 中 长 暴 有 172 个 ， 短 暴 有 8 个 。 
我 们 的 样本 是 目前 包含 缓慢 衰减 以 及 已 知 红 移 测量 的 最 大 伽 马 射线 暴 样本 ， 样 本 的 红 移 范 于 为 
(0.0368, 8.0)，7, 的 范围 为 (0. 18，549. 9) 。 样 本 中 短 伽 马 暴 的 红 移 、7T, 见 表 1。 

文 [56-57] 用 平滑 拐 折 寡 律 函 数 拟 合 这 些 具 有 平台 期 以 及 随后 正常 下 降 阶 段 的 伽 马 暴 ， 我 们 也 采 
用 该 模型 对 样本 进行 拟 合 ， 
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t =e #1 我 们 样本 中 包含 的 8 个 短 伽 马 暴 
Fy(t) = Е IG | + и =) |, (1) Table 1 8 short GRBs in our sample 

其 中 ,a 为 平台 期 下 降 指数 ， 典型 值 应 Н 0; о, СЕВ Name Т» /в 2 
描述 了 随后 的 正常 下 降 阶 段 ; T, 观察 到 平台 期 的 结束 时 СЕВ 051221A 1.40 0. 547 
间 ， 伽 马 射线 暴 静 止 系 下 的 结束 时 间 可 由 了 = 了 1 +z) 得 CRB 061201 0. 76 0. 111 
到 ; Py х 2 ”为 平台 期 末端 对 应 的 通 量 ; w 表示 光 变 曲线 СЕВ 070809 1.30 0. 220 
锐 度 参数 ， i 期 向 随后 下 降 阶 段 转变 的 锐 度 。 CRB 090510 0.30 0. 903 
我 们 采用 马尔 可 夫 链 蒙特 卡 罗 方 法 ， 以 及 卡 方 检验 获得 fan Bosh 0.18 0. 356 
We о, pre tire 吉 果 。 卡 方 检验 的 数学 公 CAR TAA aan aes 
式 为 | А CRB 150423А 0.22 1. 394 
ж = у ыбы) = Fua Ct) J А (2) СКВ 190627A 1.60 1. 942 


1 [дЕ (ti) ГА 
Жр, Е, (о) 是 观察 到 的 时 间 о 通 量 ; 6F,.(t;) 是 置信 和 度 水 平 下 的 相应 误差 。 
图 1 hi 。 通 过 拟 合 我 们 得 到 了 平台 阶段 衰减 指数  ， 以 及 缓慢 下 降 阶 段 的 误 
减 指数 a,， 平 台 期 结束 时 间 Т, 和 平台 衰减 流量 Rx ， 部 分 暴 的 拟 合 数据 见 表 2。 
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图 1 几 个 具有 平台 与 随后 正常 下 降 阶 段 的 又 射线 余辉 光 变 曲线 的 拟 合 结果 


Fig. 1 Example fitting results of X-ray afterglow light curves with platforms and subsequent normal decay phase 
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Table 2 Some fitting results and parameters of the GRBs afterglow light curve with “slow descending” 


Log Log log Log Р 
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续 表 2 
L Log log Log К 
og 2 
GRB Name Ta z 5° a =, м (Ey (Bo СЕ A 
(Т„/10°) 47) d 53 үе 535 f /dof 
10%) 10%) 10%) 
2. 3+ 0. 460+ 1. 244+ 0. 06+ 0. 3+ -1.644 -1. 63+ 
GRE doubles 2:006 0.5 0. 223 0. 183 0. 65 0. 83 0. 22 0. 22 022 
290+ 0. 986+ 1. 252+ 0. 53+ 0. 9+ 0. 11+ 0. 12+ 
СЕВ 131030A 41.1 1. 293 0 0. 285 0 036 0.18 0.31 0 0 1.04 
8. 2+ -0.243+ 1. 097+ —0. 4+ 0. 22+ -2.11+ -2. 11+ 
0 1 0. 608 0.040 1.28 0.3 0.12 0.12 с; 
71+ 0. 140+ 1. 036+ -0. 07+ 1. 09+ -0. 29+ -0. 27+ 
ЕН 112.3: 1.6% 5 0. 166 0. 049 0. 36 0. 35 0. 07 0. 07 ор 
11+ 0. 675+ 1. 314+ 1. 26+ 0. 33+ 1. 14+ 1. 12+ 
Cas ДӨ 6 2 0.077 0.811 0. 64 0.49 0.18 0.18 lade 
360+ 0. 816+ 1. 774+ 1. 18+ 0. 23+ 0. 09+ 0. 1+ 
CRB 21DolOB, 09.38. 1: 1349 10 0. 094 0. 089 0.35 0. 62 0. 03 0. 03 or 
22+ 0. 380+ 1. 658+ 1. 2+ 0. 37+ 1. 04+ 1. 04+ 
НЕ о 08 2 0. 298 0. 293 0. 46 0. 68 0. 09 0. 09 Oe 
250+ 0. 324+ 1.332+ -1. 12+ 4. 63+ 1. 0+ 1. 04+ 
СЕВ 220101A 173. 36 4.618 0 0.085 0.012 0.16 0.12 0 0 2. 48 
16+ -0.007+ 1. 639+ 0. 29+ 2. 5+ 0. 15+ 0. 1+ 
CRB 22011A.. «42-81. -4.261 2 0. 272 0. 190 0. 87 2.38 0. 12 0. 12 ә 
8. 1+ -0.326+ 1. 734+ -0. 85+ 2. 97+ -0. 37+ 0. 31+ 
GRE 0 0.8 0.147 0. 131 0.18 0.2 0.1 0.1 ИН 


注释 :单位 为 s， 取 自 Swift СЕВ 表 https ://swift. gsfc. nasa. gov/archive/ grb _table. html/. ; ba 


ПМУ 4107 erg + cm”, Ж В 


Swift GRB X https ://swift.gsfc.nasa.gov/archive/grb_table.htm//. ; “单位 为 103 s ， 用 (1) 式 拟 合 结果 ， 采 用 MCMC 算法 进 


行 拟 合 ; “单位 为 10”erg • s …， 用 (1) 式 拟 合 结果 ， 采 用 MCMC 算法 进行 拟 合 ; “单位 为 10” erg， 用 (1) 式 拟 合 结果 ， 采 
用 MCMC 算法 进行 拟 合 ; 单位 为 10” ergo 
我 们 通过 拟 合 得 到 了 F、。， 于 是 缓慢 下 降 结束 时 对 应 的 光度 为 玉 
4mDi(z) Еу ja 
а) , ) 


HEP, z: 为 红 移 ; р, 为 光度 距离 。z 和 By № Swift 伽 马 暴 数据 库 参 数 表 网 站 (https : // swift. gsfc. nasa. 
gov/archive/grb_table/) 获取 。 我 们 注意 到 在 Swift 伽 马 暴 数据 库 参 数 表 网 站 中 ， 有 时 会 给 出 两 个 红 移 
值 。 在 这 里 我 们 选取 的 红 移 值 与 文 [ 54] 一 致 。 本 文采 用 平坦 的 宇宙 学 模型 ACDM, W H, =70. 0 km · 
s' + Мре Fil Q, =0. 286 计算 光度 距离 D(z)， 


e(l +z) rf dz’ 
Diss ee Te) Н, | DT ie) FO "и 
其 中 , 0, =1-0,; c HHR, 
静止 参考 系 中 ， 观 测 到 的 瞬时 辐射 的 各 向 同性 伽 马 射线 能 量 (E') ,, ) 由 
4TD1S 
vs Tg; (5) 


计算 ， Ж, 5 为 爆发 预警 望远镜 (Swift/BAT) 通 量 ， 
fi, IRA Swift 伽 马 暴 数据 库 参 数 表 。 

通过 (5) 式 计算 Er’, ;是 不 准确 的 ， 因 为 守 宙 学 时 间 的 膨胀 ，Swift 伽 马 暴 的 观测 坐标 系 的 15 ~ 150 
keV 能 道 不 同 于 具有 红 移 的 静止 坐标 系 能 道 '*”。 因 此 ， 在 我 们 的 计算 中 应 该 考虑 所 谓 的 K-correction 
改正 。 考 虑 大 改正 效应 后 ， 修 正 后 的 伽 马 射 线 暴 各 向 同性 能 量 为 


能 量 范 围 为 (1$~ 150 keV), VA erg - cm 为 单 


202311.00002v1 
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150/1 +z 
| E®(E)dE 
15/l+z 


， (6) 


ш | “EB(E)dE 

其 中 ，BP(E) 是 伽 马 暴 的 能 谱 限 数 。 通 常 ， 件 马 暴 瞬时 辐射 的 峰值 能 量 范围 在 ~200-300 keV, НҒ 
Swift 上 爆发 预警 望远镜 探测 器 的 波段 (15~ 150 keV) 和 Swift/XRT 探测 器 的 波段 (0.2~10 keV) WEE 
范围 相对 较 窗 ， 我 们 采用 单一 的 寡 律 谱 模 型 ， 即 ФОЕ) = Е”, Ж ну 为 光谱 指数 ， 也 是 从 Swift 伽 马 暴 
数据 库 参 数 表 中 得 到 。 用 于 计算 本 文中 的 改正 的 范围 为 (1. 00, 1. 16) 。 


2 分 析 结 果 
2.1 两 类 暴 “平台 期 ”和 随后 “正常 下 降 ” 参数 的 比较 


通过 拟 合 光 变 曲 线 我 们 获得 了 两 类 暴 的 光 变 特性 参数 ， 包 括 о, a, Ao, Т,, І, M E o 927 
别 给 出 了 长 伽 马 暴 和 短 伽 马 暴 的 光 变 特性 参数 a, о, Ao ЖП Т, 的 分 布 。 


1.6 14 5 — LGRBs 


1.4 1.2 L = 


1.2 
1.0 


0.8 


Density 
Density 


0.6 
0.4 


0.2 


Density 
Density 


Aa log [T, (103s) ] 


图 2 两 类 暴光 变 特 性 参数 wj, о, AafeT, 的 分 布 ， 蓝 绿色 曲线 是 短 伽 马 暴 的 拟 合 线 ， 红 色 曲 线 是 长 伽 马 暴 的 拟 合 线 
Fig.2 Тһе distributions of two types of light curve characteristic parameters a,, a@,, Aa, and T,, where the blue-green curve 


and the red curve are the fitting lines of SGRBs and LGRBs, respectively 


通过 图 2 可 以 清晰 地 看 出 ， 对 于 我 们 样本 中 的 长 暴 和 短 暴 ， 平 台阶 段 指 数  、 正 常 衰减 阶段 指 
数 % 、 缓 慢 下 降 阶 段 和 随后 下 降 阶 段 斜 率 的 改变 量 Aa 以 及 缓慢 衰减 阶段 的 持续 时 间 Т. 都 呈现 单 峰 
分 布 。 表 3 分 别 给 出 了 两 类 暴 a, a, Aa ЯП Т, 的 分 布 范围 、 平 均值 以 及 中 值 ， 我 们 发 现 长 伽 马 暴 和 
短 伽 马 暴 的 a,, о, Aa AT, 的 平均 值 和 中 值 十 分 接近 。 
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为 了 进一步 比较 长 伽 马 暴 和 短 伽 马 暴 的 光 变 жз 长 暴 与 短 暴 平台 期 以 及 随后 下 降 各 参数 的 
特性 参数 ， 我 们 进行 了 -5 检验 ， 以 检查 这 些 SHEA, PAM FYE 
参数 的 分 布 是 否 相 同 。 表 4 列 出 了 KK-s 检验 的 Table 3 The range, median, and average values of 
结果 。 我 们 发 现 ， 对 于 о, о, Да 和 T., W 类 platform periods and subsequent decline 
暴 检验 得 到 的 值 都 大 于 10*， 由 此 我 们 可 以 判 ——_ Paneer 
定 两 类 暴 的 «|, a,, Aa 和 Т, 没有 显著 差异 。 ne NI Range Median Average 
2.2 两 类 暴 的 两 参数 关系 7.6x103s 


我 们 发 现 选取 的 样本 中 长 伽 马 暴 和 短 伽 马 暴 fa DORE igam eID O n 


《 NEF ` Sb и х. | =: А 
的 两 参数 之 间 可 能 存在 一 些 相 关 性 。 图 3 展示 SCRB #2*10 s; доруро в 68x 10° 
了 两 类 暴 的 LL MT, Ly MEY, SFY Т», Ey ico к 
和 7, 的 相关 关系 a, LGRB -0.974, 1.042 0.457 0.427 

由 图 3 可 以 看 出 ， 这 几 对 参数 的 分 布 有 -一定 SGRB 0. 122, 0.799 0.4415 0. 448 
WC, Ly -7, 的 最 佳 拟 合 关系 为 log (Ly/10" a, LGRB 0.441,4.318 1.453 1.58 
erg + s!) = (1.40 + 0.09) +(-1.44 + 0. 11) log SGRB 1.306, 2.707 1.608 1. 80 
(T,/10°s) 。 斯 皮尔 曼 (Spearman ) 相关 系数 为 -0.72， Aa LGRB 0.92, 1.91 1.02 1.16 
置信 概率 p < 10”。 我 们 可 以 看 出 ， 两 类 暴 缓慢 SCRB 0.0, 4.83 1. 365 1.35 
衰减 阶段 的 持续 时 间 Т, 与 缓慢 下 降 阶 段 结束 时 
сак жа 两 类 暴 缓慢 下 降 以 及 随后 下 降 各 参数 的 
БН KAA log( Ly/10* erg + в) = (0.48 + K-S 检验 也 值 以 及 p 值 
0.06)+(1.16+0.06)log( Е 10е) 。 斯 皮尔 Table 4 K-S test results of the two types of GRBs with 
曼 相关 系数 为 0.739， 置 信 概 率 p=2. 408x10”, slow decline and subsequent decline 
因此 ， 两 类 暴 缓慢 下 降 阶段 结束 时 刻 对 应 的 光度 K-S Test D value E 

上 p Als -E yi АЕ 
Ly 3 4 Fl [F] PE RE Bt iso = МИТЕ FE IEMA AER o Т, 0.2326 0. 725 
除 此 以 外 ， 我 们 发 现 两 类 暴 的 $, Ey ТЛЕ Р КЕРЕ нет 
正 相 关 关 系 。 其 中 ，S-7w 的 最 佳 拟 合 结果 为 log | | 
а, 0. 366 3 0. 202 


(85/107 erg • стг) = (0. 45 + 0. 10) +(0. 52 + 0. 05) 
log( To/s) ， 斯 皮尔 曼 相 关系 数 为 0. 57， 置 信 概 个 0470 о% o 
R p<10*, Е, „ТН 0 log( E,/ 

10° erg) = (-1.14+0. 15) +(0. 78 +0. 15)Ty)/s), BTR A SH KARA 0.47, BAE p<10%. Н sS- 
ТЖЕ, -7 的 关系 可 以 看 出 ， 观 测 到 的 伽 马 射线 暴 的 持续 时 间 越 长 ， 能 量 / 能 流 越 高 。 

此 外 ， 在 图 3 中 我 们 用 黑 圆 点 标记 长 伽 马 暴 ， 用 红 圆 点 标记 短 伽 马 暴 ， 以 便 观察 这 两 类 伽 马 暴 是 
否 存 在 差异 。 从 图 3 KNITS, KMS LMT, LAE, io SAT), Ey „ЛП Too 
的 相关 关系 是 一 致 的 ， 基 本 分 布 在 最 佳 拟 合 结果 附近 。 由 于 两 类 暴 的 ЯП Т, LA E io ЗАТ, 
Esco FU Tom 的 相关 关系 没有 呈现 明显 的 差异 ， 因 此 我 们 推测 ， 这 两 类 伽 马 暴 在 本 质 上 可 能 是 相同 的 一 
类 暴 ， 因 而 它们 的 喷 流 结构 也 可 能 相同 。 


3 讨 论 


由 于 伽 马 暴 的 分 类 和 物理 起 源 等 问题 还 不 清楚 ,本文 主要 从 伽 马 暴 X 射线 余辉 的 时 间 特 性 出 发 ， 
通过 使 用 平滑 拐 折 寡 律 函数 拟 合 具有 缓慢 下 降 以 及 随后 正常 下 降 的 两 类 伽 马 暴 的 X 射线 余辉 光 变 曲 
线 ， 得 出 并 比较 两 类 暴 的 X 射线 余辉 时 间 特 性 参数 ， 从 而 推断 这 两 类 伽 马 暴 喷 流 结构 的 物理 起 源 是 
和 否 一 致 。 我 们 选择 180 个 包含 已 知 红 移 的 伽 马 暴 样 本 ， 其 中 长 伽 马 暴 172 个 ， 短 伽 马 暴 8 个 。 我 们 使 
用 马尔 可 夫 链 蒙特 卡 罗 方 法 以 及 卡 方 检验 对 这 些 包含 平台 的 余辉 数据 进行 拟 合 ， 从 而 获得 了 光 变 特性 
BH о, а, ЯП Т, 并 通过 计算 得 出 L M E ioo 
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一 Best fit line 一 Best fit line 
6 в Long GRBs in our sample 6 в Long GRBs іп our sample 


ж Short GRBs in our sample ж Short GRBs in our sample 


log (L, / 1047 erg . 571) 
log (L, / 1047 erg . s!) 


log (Т, / 10° s) log (Ey, iso / 10° erg) 
(a) (b) 
— Best fit line — Best fit line 
3.0 в Long GRBs іп our sample = 27 f Long GRBs іп our sample 
ж Short GRBs in our sample ® Short GRBs in our sample I т ж i 
25 x ~ I ғ if pa т 
р ъй 
5 2.0 S ob 
5 15 = 
ра зт 
10 外 
2 > 
® 05 = -2p 
0.0 -3 上 ғ? I 
| I 
-0.5 —1- L ee ee НЫНЕ ee See eee 
-0.5 0.0 0.5 1.0 1.5 2.0 25 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 
log (То / 8) log (Тоо / $) 
(с) (а) 


图 3 两 类 暴 的 两 参数 关系 的 对 数 分 布 图 。(a)L-7T,; (b)Ly-E iw; (c)S-7T%; (ADE, 
红 圆 点 是 短 伽 马 暴 ， 蓝 色 实 线 为 最 佳 拟 合 线 
Fig. З Two parameters relationship diagram for two types of GRBs. (a)Ly-T,; (b)L,-E,, 


„То. BARKER, 


„iso 


3 (c)S-Tyy; (d)E, sTo The 
black and red dots are the LGRBs and SGRBs, respectively, and the solid blue lines are the best fit line 


iso 


我 们 首先 比较 了 长 短 暴 X 射线 余辉 光 变 曲线 的 参数 a, а, ЛП Т, 是 否 存在 差异 ， 通 过 对 两 类 暴 的 
均值 、 中 值 以 及 K-S 检验 后 发 现 ， 两 类 暴 的 光 变 特性 参数 分 布 没 有 显著 差异 。 

我 们 也 对 比 了 静止 坐标 系 下 180 个 伽 马 暴 的 缓慢 衰减 结束 时 间 7,、 平 台 结束 时 对 应 的 光度 工 ,、 
各 向 同性 能 量 瑟 ,的 分 布 情况 。 我 们 的 研究 发 现 ， 缓 慢 下 降 结 束 时 间 7,、 平 台 结 束 时 对 应 的 光度 工 、 
各 向 同性 能 量 5, ,的 分 布 都 接近 一 个 单 峰 分 布 。 这 些 参 数 的 分 布 范 围 与 文 [| 50] 和 文 [54] 样本 的 参数 
结果 非常 接近 ， 说 明 我 们 的 分 析 结 果 和 前 人 的 结果 一 致 。 

文 [47] 分 析 了 77 个 伽 马 暴 ， 发 现在 静止 坐标 系 中 缓慢 衰减 阶段 的 结束 时 间 (7,,) 与 对 应 的 X 射 
线 光 度 ( 工 , ) 之 间 存 在 负 相 关 关系 ， 即 оет, O”, 我们 的 研究 发 现 ， 当 使 用 一 个 更 大 的 样本 重 
新 分 析 这 两 个 量 的 相关 性 时 ， 两 参数 相关 性 也 存在 ， 斯 皮尔 曼 相关 系数 p 为 -0.72， 和 置信 概率 p< 
10“， 可 以 看 出 该 相关 性 是 很 强 的 。 随 后 ， 文 [49] 通 过 已 知 红 移 的 62 个 长 Swift WAHRER, H Ly ,- 
7, ,的 相关 性 扩展 到 瞬时 辐射 之 间 ， 并 发 现 Ly ,与 瞬时 辐射 平均 光度 Lo 5 存在 相关 性 ， 即 = Ls 
[50] 进一步 发 现 了 ,， 太 ,和 了 瞬时 辐射 的 各 向 同性 能 量 (已 ,,) 之 间 存 在 更 紧密 的 三 参数 相关 性 ( 简称 
为 L-T-E 关系 ) 。 文 [51] 建 立 了 瞬时 辐射 中 了 ,， 太 ,和 峰值 光度 ( 户 ) 之 间 的 另 一 个 三 参数 相关 性 。 
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文 [53] 分 析 了 50 个 在 光学 余辉 光 变 曲线 中 具有 缓慢 衰减 阶段 的 伽 马 暴 ， 发 现 样 本 具有 相似 的 三 参数 
L-T-E 相关 性 。 后 来 ， 文 [54] 分 析 了 174 个 具有 缓慢 衰减 阶段 和 已 知 红 移 的 Swift 伽 马 暴 大 样本 ， 并 
测试 了 L-T-E 的 相关 性 。 

从 理论 上 看 ， 目 前 文献 中 提出 了 大 量 的 模型 试图 解释 X 射线 平台 。 文 [60] 提 出 了 伽 马 暴 的 毫秒 
原 磁 星 模型 ， 在 该 磁 星 模型 中 ， 偶 极 子 辐射 释放 的 能 量 预 计 通 过 磁 星 风 传递 到 周围 环境 。 这 种 风 是 大 
质量 恒星 核心 十 缩 产生 超新星 激 波 后 ， 由 中 微 子 加 热 形 成 的 。 通 过 与 前 身 星 的 相互 作用 ,流出 物 被 准 
直 成 双 极 喷 流 。 随 着 磁 星 冷 却 ， 变 成 极端 相对 论 ， 坡 印 享 通 量 在 时 间 尺 度 上 占 主导 地 位 ， 时 间 尺 度 与 
喷 流 清除 恒星 空 腔 所 需 的 时 间 相 当 。 因 此 ， 在 这 个 模型 中 ， 磁 耗 散 和 激 波 解释 了 瞬时 辐射 和 随后 的 陡 
峭 衰 减 阶段 。 晚 期 的 炊 发 或 余辉 辐射 ， 例 如 X 射线 平台 ， 被 认为 是 由 持续 的 磁场 或 剩余 的 旋转 能 量 
提供 动力 。 事 实 上 , 文 [61] 用 这 个 模型 解释 了 le 不 确定 性 内 的 光度 和 平台 期 持续 时 间 之 间 的 相关 
性 ， 其 中 磁场 和 快速 旋转 中 子 星 的 自转 周期 是 典型 的 。 随 后 ， 文 [62] 的 研究 表明 ， 只 有 在 具有 高 磁 
场 强 度 的 超 磁 星 中 ， 磁 星 模型 才能 在 件 马 射线 暴 平 台 辐 射 内 一 致 。 

文 163] 首 次 将 自转 磁 星 的 非 理 想 模 型 拟 合 到 余辉 数据 中 ,统计 样 本 为 40 个 具有 明确 平台 的 长 伽 
马 射线 暴 和 13 个 短 伽 马 射 线 暴 ， 包 括 具 有 扩展 辐射 的 短 伽 马 射线 暴 ， 得 出 的 结论 是 短 伽 马 射线 暴 和 
长 伽 马 射线 暴 都 可 以 用 磁 星 模型 解释 。 在 本 文中 长 暴 和 短 暴 的 a 和 a, 分 布 都 是 单 峰 分 布 ， 而 且 两 类 
ЖЕН] a, Ma, 分 布 没 有 明显 的 差异 。 同 样 的， 长 又 和 短 暴 的 7, 和 工分 布 都 是 单 峰 分 布 且 无 明显 差 
异 。 这 些 结果 似乎 表明 ， 两 类 暴 是 同一 起 源 且 都 可 以 用 磁 星 模 型 解释 。 

文 [64] 显 示 了 在 一 些 短 伽 马 射 线 暴 中 看 到 的 扩展 辐射 与 平台 期 之 间 的 联系 。 他 们 假设 一 个 来 自 
中 心 致 密 物 体 (如 磁 星 ) 的 持续 能 量 注 和 人 ， 作 为 两 个 辐射 阶段 的 起 源 。 在 一 些 长 伽 马 射线 暴 中 也 报道 
了 这 种 持续 10° ~10°s 的 扩展 辐射 。 文 [65] 认 为 ， 在 少数 情况 下 ,平台 期 之 后 会 出 现 急剧 衰减 ， 之 后 
可 能 出 现 第 2 个 平台 期 。 他 们 将 此 解释 为 中 央 引 警 的 结果 ， 该 引擎 由 一 个 快速 旋转 的 磁 星 组 成 ， 该 磁 
РАА УТЕ Al 

与 之 前 的 工作 相反 ， 文 [66] 认 为 ， 如 果 新 生 磁 星 的 磁场 处 于 瞬 态 的 快速 衰减 阶段 ， 则 可 能 无 法 
实现 清晰 的 平台 阶段 。 由 于 这 个 阶段 自转 减 慢 ， 功 率 可 能 会 下 降 得 较 快 ， 因 此 导致 缺乏 平稳 阶段 。 基 
于 这 种 想法 ， 他 们 分 析 了 6 个 伽 马 射线 暴 的 X 射线 光 变 曲线 没有 平台 ， 并 得 出 结论 ， 这 些 伽 马 射线 
暴 宿主 可 能 是 毫秒 磁 星 。 

文 167-68] 提 出 了 一 种 不 同类 型 的 模型 ， 其 中 平台 期 、 兆 发 和 可 能 的 陡坡 都 是 由 反 向 激 波 主导 解 
释 的 '21  。 这 种 反 向 激 波 不 同 于 “火球 ”模型 中 预测 的 经 典 反 向 激 波 ， 因 为 它 被 认为 是 通过 极端 相对 
论 性 的 、 高 磁化 的 脉冲 星 状 风 传 播 的 ， 这 种 风 由 持久 的 中 央 发 动机 产生 。 该 想法 虽然 有 趣 ， 但 到 目前 
为 止 ， 只 有 极 少 数 观 测 到 的 伽 马 射线 暴 似乎 符合 这 个 模型 的 预测 。 此 外 ， 这 些 伽 马 暴 在 其 X 射线 余 
辉 中 有 多 个 炮 发 ， 因 此 很 难得 出 确切 的 结论 。 在 经 典 “ 火 球 ” 模 型 的 范围 内 “7 ， 从 滑行 到 相似 阶 
段 的 过 渡 过 程 中 ， 反 向 激 波 的 贡献 是 预期 的 。 事 实 上 ， 在 文献 中 我 们 确实 看 到 了 一 些 伽 马 射线 暴 的 光 
变 曲线 可 能 是 来 自 反 向 激 波 的 证 据 ( 如 GRB 110213A)。 

文 [72-73] 提 出 的 另 一 个 模型 是 从 离 轴 角 度 看 结构 喷 流 。 这 些 工作 的 目的 是 从 极 低 光 度 GRB/GW 
170817A 中 检测 离 轴 喷 流 。 文 [72] 通 过 高 纬度 辐射 解释 了 平台 期 ， 文 [73 ] 认为 在 近 核 心 视线 上 存在 
结构 化 喷 流 。 因 此 这 些 平台 将 会 被 大 量 的 观察 者 看 到 ， 并 持续 10° ~ 10s。 它 自然 地 再 现 了 观测 到 的 
时 间 尺 度 和 光度 分 布 ， 以 及 平台 持续 时 间 、 平 台 光度 和 瞬时 y 射线 能 量 之 间 的 相关 性 。 在 本 文中 两 
类 暴 也 存在 L-7, ВІ, -Е, ;相关 性 ， 这 可 能 表明 两 类 暴 是 同一 起 源 ， 因 而 可 以 推断 它们 的 物理 机 制 
应 该 相同 。 

由 本 文 得 出 两 类 暴 的 光 变 特性 参数 具有 相似 的 分 布 ， 且 文 [63] 得 出 的 结论 是 具有 缓慢 下 降 阶 段 
的 两 类 暴 可 以 在 磁 星 模型 中 解释 ， 文 [73] 提 出 的 模型 自然 地 解释 了 平台 的 持续 时 间 和 光度 ， 以 及 瞬 
时 y 射线 能 量 之 间 的 相互 相关 性 。 这 些 证 据 表明 ， 长 暴 和 短 暴 可 能 是 同一 起 源 。 

因此 ， 如 果 我 们 想 要 了 解 更 多 长 暴 和 短 暴 的 性 质 ， 还 需要 更 大 的 样本 进行 深入 研究 。 
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AS SCE AES АЛ FEE PROT DA 218 BEL, ЖЕ Г БЕЛИЛ BY DS eS EE АЧТА AE 
曲线 分 析 。 通 过 对 X 射线 余辉 光 变 特性 分 析 发 现 ， 两 类 暴 的 光 变 曲线 中 都 含有 缓慢 下 降 以 及 随后 下 
降 的 成 分 。 我 们 对 比分 析 了 两 类 暴 的 光 变 特性 参数 得 到 以 下 有 趣 的 结论 : 

(1) 我 们 选择 的 两 类 暴 的 Х 射线 余辉 光 变 特性 参数 a, a, Aa ЯП Т, 的 分 布 都 是 单 峰 的 。 

(2) 我 们 得 出 了 两 类 暴光 变 特性 参数 的 范围 、 中 值 和 均值 : af(-0.974，1. 042) 0. 457, 0.427; 
ai(0. 122, 0.799), 0.441 5, 0.448; ay(0.441, 4.318), 1.453, 1.58; af(1.306, 2.707), 1.608, 
1.80; 7'(0. 85x 10°s, 3. 68 х 10's) , 4.91 10°s, 8.68X10°s; 以 及 两 类 暴光 变 特 性 参数 的 K-S 检验 p 
ЇН, a, 50.977, о, 40.202, Т, 为 0.725， 其 值 均 大 于 10“， 这 些 值 表 明 两 类 暴 的 光 变 特性 参数 分 
布 没 有 显著 差异 。 

(3) 我 们 还 发 现 X 射线 缓慢 下 降 阶 段 的 拐 折 时 间 从 几 十 秒 到 10°s 不 等 ， 典 型 值 约 为 10's， 对 应 的 
FEN (Ly ) 主要 在 10“ erg +5 A] 10% erg .ss 之 间 ， 一 般 比 光学 余辉 平台 对 应 的 光度 要 高 两 到 三 个 数量 
级 ， 这 支持 了 文 [53] 的 观点 。 

(4) 同时 我 们 还 发 现 选取 的 两 类 暴 样本 中 存在 一 些 有 趣 的 两 参数 相关 性 如 L, -T,(Dainotti KA ) : 
log(ZL AL10”erg.s)=(1.40+0.09)+(-1.44+0.11)log(7.Z102s) ， 其 斯 皮尔 曼 相 关系 数 为 -0.72， 置 
信和 概率 p=1.42x10™; І.-Е: 1ов(1,/10”егв+ 5 ')= (0. 48 + 0. 06+( 1. 16 + 0. 06) log( E,/10” erg) , 
斯 皮尔 曼 相 关系 数 p 为 0.739， 置 信 概 率 为 p=2. 408 x 10°; S-Ty: log(S/10erg + em”) = (0. 45 + 
0. 10)+(0.52+0.05)log(7TuXs) ， 斯 皮尔 曼 相 关系 数 p 40.57, В p=5.90x 10°"; Е, 
Т»: log( E,/10° erg) = (—1. 14+ 0. 15) +(0. 78 + 0. 15) log( Ты/з), Spearman 相关 系数 p 为 0.47， 置 信 
概率 为 p<10“， 且 两 类 暴 的 两 参数 相关 关系 没有 显著 差异 。 

以 上 的 研究 结果 表明 ， 两 类 暴 的 X 射线 余辉 光 变 特性 在 参数 分 布 和 两 参数 相关 关系 上 具有 相似 
的 特征 ， 这 也 为 两 类 伽 马 暴 可 能 是 同一 物理 起 源 提 供 支 持 。 
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Abstract: The classification of gamma-ray bursts (GRBs) is an unsolved hot issue in the field of GRBs. 
To solve this problem, we collected a samples of X-ray afterglow curves including 180 known redshift GRBs 
observed by Swift with slow attenuation and subsequent attenuation characteristics, in which there are 172 
Long GRBs and 8 Short GRBs. We adopt the smooth broken power law function and the Markov chain Monte 
Carlo (MCMC) method to fit these X-ray afterglow light curves. After obtaining these light curve characteristic 
parameters a, (slow decay index), œ, (normal attenuation index), and T, (slow attenuation end time in rest 
frame), we analyze the characteristic parameters of the two types of bursts to check whether there are 
differences in the distribution of these parameters and whether the two types of storms follow the same 
correlation. Our research findings are as follows: (1) The distribution range, median value, and average value 
of the characteristic parameters (a@,, @,, T,) of two types of bursts are very close, and there is no significant 
difference. The subsequent K-S test further indicates that there is no significant difference in the distribution of 
parameters between the two types of bursts. (2) The two types of bursts we selected follow the following 
S-T,,. Ё 


statistical results indicate that the optical variability characteristics of the two types of bursts are similar. 


correlation; L,-T,, LE 5... уло Тоо, and they аге not distinguishable on a scatter chart. These 


Therefore, we infer that the two types of GRBs have the same physical origin, and their jet structures may also 
be the same. 
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